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1. Introduction 

Lattice QCD with chiral fermions, although computationally expensive, is the best formulation 
of QCD on the lattice. There are two chiral formulations: a) Domain wall fermions [jlj, ||] and b) 
overlap fermions [§], which are closely related [Q |5j|. 

In this work we focus on the use of domain wall type fermions for lattice QCD simulations. 
Recent dynamical simulations with such fermions revealed larger chiral symmetry breaking than 
expected [B]. It was shown that the effect of gauge action improvement is effectively canceled 
by the dynamical domain wall fermion. Since there is no such problem with dynamical overlap 
fermions one wonders if there is something wrong with domain wall fermion formulation. 

But as mentioned above these formulations are closely related and therefore if anything is to 
be blamed it is the implementation of domain wall fermions. It is the purpose of this work to show 
that the state-of-the-art of the dynamical domain wall implementation is not well suited for the 
state-of-the-art simulation algorithms. 



2. Notations 



In this work we will use the truncated overlap fermions [f7|] which have better chiral properties 
than the standard domain wall fermions The corresponding 5-dimensional operator is given by 
the N5 x N5 blocked matrix: 



( D w -l 
(D W + 1)P_ 

\-m{D w + l)P+ 



(D w +l)P + 
D w -1 



-m 



(D w + l)P\ 



(D w + l)P+ 
(Dw + 1)P- D w -l j 



where the blocks are matrices defined on the 4-dimensional lattices using the negative mass Wilson- 
Dirac operator Dw, N5 is the number of time slices along the fifth Euclidean dimension and P± are 
chiral projection operators. Note that standard domain wall fermions use off-diagonal blocks which 
omit D w . 

Let ^#1 be the same matrix as above but with the bare quark mass m= 1. Then it can be shown 
thatfi|]: 

(2.1) 



det^f 1 ^ = detD (A ' 5) 



where 



fN 5 



2 1 2 ,J l + 7*& 
with T the transfer matrix along the fifth dimension given by: 

^ l + H w 



(2.2) 



w 



H w = YsDw 



Note also that in the large limit approaches the Neuberger overlap operator 

^ 1 +ra 1 —m . . . 
D = —— 1 + — — sign (Hw ) 
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3. The problem 

The standard way one deals with the fermion determinant is expressing it as a Gaussian integral 
over pseudofermion fields: 



|detD^| 2 = \fetJt{ l J(\ 2 = J [d&d$\e 



if 



The resulting effective fermion action is given by: 

Spf = 



"^lOll 2 



From equation 2.1 we expect the action to have one contribution from the operator and 
some extra terms: 



Spf 



Xi 



N 5 

l 2 + Ll|crf 



where Xi pseudofermion fields and C, are four dimensional matrices which in general may 
depend on gauge fields. This form of the action will be explicitly calculated below in this paper. 

We note the extra terms in the action and we ask whether they would contribute in the gen- 
eration of gauge fields. As long as the simulation algorithm averages over a large ensemble of 
pseudofermion fields Xi these terms will cancel to give the correct determinant. However, fermion 
algorithms typically construct molecular dynamics trajectories which keep the pseudofermion field 
fixed. This may cause the exploration of gauge field configurations which "feel" the extra dimen- 
sion through the extra terms. 

These extra pseudofermion terms may be regarded as "artefacts" of the algorithm which con- 
tribute via the renormalised coupling: 

j3^ Cl j3+Aj3 

As observed by |Q] the renormalisation is such that it drives the gauge field toward the Aoki parity 
broken phase which in turn causes the breaking of the chiral symmetry. 

We believe that this effect should go away if the extra terms are not present in the action. A 
direct evidence for this is not provided here, but simulations with overlap fermions clearly show 
that such renormalisation effects are absent. 



4. The solution 

First, let us calculate 5pp. Using algebraic manipulations as in [||] and the results of the Ap- 
pendix we get: 

s PF = H^-^oH 2 = \\^- x 3r x p T ^\\ 2 

where 2T is given by: 

/ P + -mP -T \ 



\-T(P--mP+) 



-T 
1 / 
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and P is the permutation matrix: 

(P+P- 
P + '■ 
'-.P. 
\P- P+J 
It is straightforward to show that (see Appendix): 

/[#)]-' 

S (2) 



\ 1 J 

where flW = T N ^ i+1 {(P_ - mP + )[D^}- 1 + 7s = 2, . . . ,N 5 . This way we get: 

A'5 



(4.1) 



5 PF = 1 1 1 1 2 + £ | |fl(04>! + | 2 



i=2 

where 4>i = P + <3>i + P <3>/v 5 and 4>; = P + <E>, +P <!>, i , i = 2, . . . ,Ns. It is easy to see that the second 
term of the right hand side constitutes the bulk of the coupling constant renormalisation. 

This result hints to the following solution of the problem: express the fermion determinant in 
such a way that there are no additional terms in the effective fermion action. An obvious solution 
is to define: 

5 PF = ||[D^]->|| 2 



where <j> is now the usual pseudofermion field defined on the four dimensional lattice. The difficulty 



with this form is that the differentiation of ( |2.2[ ) w.r.t. gauge field may yield numerically unstable 
expressions and ill-conditioned matrices. 

A stable implementation may be defined as in the following: let E\ be denned as the blocked 
unit direction along the fifth dimension, i.e., 

£l = (l,0,...,0) r 



Using ( |4.1| ) we can write for the fermion determinant: 

|detD (iV5) | 2 = Idetef^f 1 ^! 2 = (detefP^f^Peil 2 

This expression can be used as a starting point to formulate a simulation algorithm in terms of 
numerically stable derivatives of ^-matrices. Hence, a variation of can be computed using 
the variations of and . 

Note that the equation \.\ can be used to compute the inverse of D^ Ni \ In this case the linear 
system to be solved is: 

D {Ns) x = b (4.2) 
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Let y = (x, yi ■ ■ ■ ,yN 5 ) and z = (b,0...,0) be 5-dimensional vectors. Then, from eq. 4.1 (see 
Appendix) we can write: 

3Ty = ZT x z 

Using eq. [TT] we get: 

In this way, one can use the usual solver for the 5-dimensional matrix and get the solution x of the 
4-dimensional system 1.2 form the first block-component of y. 



5. Conclusion 

In this paper we have given an explicit calculation of the effective pseudofermion action that 
is used in dynamical domain wall simulations. We note that the extra pseudofermion terms may 
cause the current simulation algorithms to produce gauge fields in the region of the Aoki phase. 

To cure this phenomenon we have proposed a simple solution which can be easily implemented 
using the algebraic relations between overlap and domain wall fermions. 
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6. Appendix: explicit calculation of Spp 

Multiplying ^# from the right by P we obtain: 



75 



f(H w -l)(P + -mP-) H w + l 

H w -l'- 



'■■ H w + l 
H W -1J 



\(H w + l)(P_-mP + ) 
Further, multiplying this result from the left by the inverse of the diagonal matrix: 

H 5 = y 5 dmg(H w -l,... : H w -l) 

we get: 

3? = H 5 l JKP 



(6.1) 



hetXjj,i,j = 1, . . . ,Ns be the four dimensional blocks of the inverse of 2F . Then a straightfor- 
ward calculation gives: 



X, 



lU (P+ - mP ) - T N * (p_ - m p 
for j = 1, . . . ,iVs and 



i-j 



-75- 



-l 



for i = 2,... ,Ns and 



= (P_ - mP + )X xj + £ S^jT* 

for i = 2,...,N 5 , j = 2,...,N 5 . 



N 5 -i 

I 

k=0 
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